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C ■ ABSTRACT 

^^ ■ In previous work, we have shown that recent updated standard solar models cannot reproduce 

I , the radial profile of the sound speed at the base of the convective zone and fail to predict the 

pq ■ photospheric lithium abundance. In parallel, helioseismology has shown that the transition 
from differential rotation in the convective zone to almost uniform rotation in the radiative 



solar interior occurs in a shallow layer called the tachocline. This layer is presumably the 

seat of a large scale circulation and of turbulent motions. Here, we introduce a macroscopic 

Xq I transport term in the structure equations, which is based on a hydro dynamical description of 

ff^ . the tachocline proposed by Spiegel & Zahn, and we calculate the mixing induced within this 

^O \ layer. We discuss the influence of different parameters that represent the tachocline thickness, 

the Brunt- Vaissala frequency at the base of the convective zone, and the time dependence of 

Qs \ this mixing process along the Sun's evolution. We show that the introduction of such a process 

inhibits the microscopic diffusion by about 25 %. Starting from models including a pre-main 

sequence evolution, we obtain: 

a) a good agreement with observed photospheric chemical abundance of light elements such 

g : as ^He, ^He, ^Li and ^Be, 

^1 b) a smooth composition gradient at the base of the convective zone, and 

C^ ■ c) a significant improvement of the sound speed square difference between the seismic Sun and 

the models in this transition region, when we allow the photospheric heavy element abundance 
to adjust, within the observational incertitude, due to the action of this mixing process. 
The impact on neutrino predictions is also discussed. 
C^ I Subject headings: Sun: interior. Sun: tachocline. Sun: oscillations, abundances, lithium, 

neutrinos 

1. INTRODUCTION gen burning phase (Tomczyk et al. 1995, Kosovishev 

The validation of solar structure is crucial to check ^^ al. 1997, Turck-Chieze et al. 1997, Christensen- 

the hypotheses used in stellar evolution modelling. ^ ^' \ 

It allows also the helium content at the birth of the ^^ V^^^\\^\ the development of neutrino detectors 

Sun to be determined, which is key to galactic enrich- ^^ P^^icle physics has the_ objective of determimng 

ment and for predicting the present neutrino emission n i r- ^ ' r 

fluxes. With the support of hehoseismology, we are ^^^ P^^^^' °""« ^S^^"^' ^ fPf ^^^ \''}\'^ ^ ^°f "" °f 

now able to constrain the sophisticated calculations electronic neutrinos v which could be transformed 

performed to describe the present structure of the mto other flavors before reachmg the detectors on 

Sun. Very recently, a tremendous effort using ground ^arth (Hata & Langacker 1997, Bahcah, Krastev & 

based networks and space experiments has been made Smirnov 1998) . The solar conditions are really mac- 

to improve the accuracy of the seismic tool at a level f^^^^le in experiments on Earth, so one must pre- 

where differences with the theoretical estimates are ^^^^ precisely the emitted neutrino fluxes m order to 

significant. Consequently, it has been demonstrated P'jJ'^'' '"^ °''' knowledge of these particles and the 
that this discipline has the potential to check classi- 

cal stellar structure for solar-type stars in the hydro- ^* ^^^^ S^^^^*^' hehoseismology shows that the stan- 



dard solar framework is correct (Brun, Turck-Chieze 
&: Morel 1998; hereafter paper I). It is nevertheless 
extremely clear that standard models have always 
failed to answer some well known questions or anoma- 
lies due to its simplified assumptions (Turck-Chieze 
et al. 1993). These open questions concern the ef- 
fect of the rotation and magnetic field, the history of 
the angular momentum, the real way the convection 
acts, the unexplained photospheric abundances as the 
photospheric lithium depletion and the presence of 
mixing at different locations inside the stars. Recent 
results obtained with the SOHO satellite (Frohlich 
et al., 1997, Kosovishev et al. 1997, Gabriel et al. 
1997) and ground networks have demonstrated that 
there exists significant differences in the square of the 
sound speed between the Sun and the models (Koso- 
vishev et al. 1997, Turck-Chieze et al. 1997). Some 
of them can be attributed to uncertainties in the mi- 
croscopic description of the Sun, while some others 
must be the sign of missing processes in the standard 
solar framework. 

The idea of this paper is to introduce a new term 
in the stellar structure equations in order to progress 
beyond the standard framework. This new term de- 
scribes some physical processes which are presumably 
present in our star as suggested by the differential ro- 
tation profile. In this study, we hope to treat more 
correctly the transition layer between the convective 
and radiative zones, with the help of crucial measured 
helioseismic variables. The main reasons for our in- 
terest in this layer are: 

• The change of rotation rate (solid =^ differen- 
tial) measured by helioseismology (Thompson 
et al. 1996, Kosovichev et al. 1997) devel- 
ops probably some turbulent mixing in a shal- 
low layer (tachocline) which could be where 
the magnetic field is anchored and amplified 
by dynamo process (Schiissler 1987, Choudhuri, 
Schiissler &: Dikpati 1997), 

• The fact that the standard solar model frame- 
work including microscopic diffusion and ele- 
ment settling has failed to treat this transition 
correctly (paper I). Such a model leads to sharp 
element profiles which create a large bump in 
the sound speed profile below the convective 
zone and cannot reproduce the observed lithium 
depletion, 

• The turbulent convective zone surely interacts 
with the stably stratified radiative zone (e.g the 
Brunt- Vaisala frequency N"^ > 0), developing 
turbulence, modifying the thermal structure of 
the star (convective penetration) (Zahn 1991), 
or extending the mixing zone (overshooting) 



(Roxburgh 1997 and Zahn 1998) or generating 
internal waves (Press 1981, Schatzman 1993). 

Following the description of Spiegel & Zahn (1992) 
(hereafter; SZ92) of the motions in this so-called 
tachocline and that of Chaboyer & Zahn (1992) for 
the chemical evolution, we introduce in the equations 
of the stellar structure an effective diffusion coeffi- 
cient and build several solar models, which are com- 
pared to the helioseismic and surface abundance ob- 
servations. 

We first present in section 2 the physical model of 
the solar tachocline we use and deduce the vertical 
macroscopic diffusion coefficient to be introduced in 
the time evolution equation of the chemical abun- 
dances. Then in section 3 we present solar models 
with a mixing occuring in the tachocline, its impact 
on photospheric abundances, especially for lithium 
and beryllium and discuss the effect of the time evo- 
lution of such a mixing from the Sun early phases 
until the present days. In section 4 we comment our 
results on the changes observed in the sound speed 
profile, on the role of the solar composition and the 
recent nuclear reacton rates and on neutrino predic- 
tions. Finally we conclude in section 5. 

2. THE SOLAR TACHOCLINE 

The Sun is generally assumed to be in hydrostatic 
and thermal equilibrium, neglecting the effects of ro- 
tation and magnetic field. However, models built 
with these simplifying assumptions do not totally 
agree with the helioseismic data, in particular with 
those obtained by the satellite SOHO. It appears that 
macroscopic mixing processes must be taken into ac- 
count not only in the convective zone, but also in 
the radiative interior. Macroscopic mixing may be 
introduced in the solar models by adding an effec- 
tive diffusivity Dt in the equation governing the time 
evolution of the chemical species Xi. The equation 
becomes: 



dt 



dm 



+ nuclear terms, (1) 



where the velocity Vj of species i with respect to the 
center of mass is: 



V^ = -4^Tpr\D, + DT) 



dlnXi 
dm 



+ Vi. 



(2) 



The velocity Vi is the sum of a term which de- 
pends on the concentration gradient and another 
which does not, Vi (Burgers 1969, Proffit & Michaud 
1991). In this paper, we add a macroscopic diffu- 
sion term Dt to the microscopic diffusion term Di. 
This will describe the mixing occurring in the shear 



layer connecting the differential rotation in the con- 
vection zone to the solid rotation in the radiative inte- 
rior. As presented in SZ92, a possible physical inter- 
pretation of this tachocline is obtained by invoking 
anisotropic turbulence with much stronger viscous 
transport in the horizontal than in the vertical di- 
rection (i.e uh » ^y)- Such turbulence reduces the 
differential rotation and therefore inhibits the spread 
of the layer deep inside the radiative zone, leading to a 
stationary state. For a justification of this anisotropy, 
we refer to Michaud & Zahn (1998). 

2.1. The hydrodynamical description 

The conservation of mass, momentum and entropy 
is given by: 

dp 



dt 



+ V.{pV) = 



(3) 



dV 



/9 I — + {V.V)V + 2nAV + nAr\ = -VP 



-/jV$ + v. II r II (4) 



PT[J^+V.V 



s = v.(xvr) 



(5) 



where p is the density, P the pressure, T the tem- 
perature, V = {u,v,rQsm9) is the local velocity in 
the rotating frame fi. In this expression, the differ- 
ential rotation Q appears explicitly, S is the specific 
entropy, || r || the viscous stress tensor and $ the 
gravitational potential. In order to ease the solution 
of this system of equations, some simplifying assump- 
tions were made by SZ92, mainly: 

• velocities are small compared to the sound 
speed: dp/dt is negligible (anelastic approxi- 
mation) , 

• advection terms and viscous forces are small 
compared to the Coriolis acceleration (i.e small 
Rossby and Ekman numbers), 

• the tachocline is thin compared to the pres- 
sure scale height and the turbulent horizontal 
viscosity is stronger than the vertical one (i.e 
i^H » i^v)- 

After separating each variable into its mean value 
on the sphere plus a perturbation, e.g. T{r^ t) + 
T(r,9,t), the linearized form of the equations in a 
stationary state is given by SZ92. Spiegel and Zahn 
define a streamfunction ^ for the meridional flows by: 



r pu 



dx 



rpv sm ( 



dr 



(6) 



where x = cos 6, and project the variables of the 
equations describing the tachocline on horizontal 
eigenfunctions -Fj: 



{P,T,u) 
xil 



^{Pi,fi,Ui)Fi{x 

i 

J2^i f Fi{x)dx 



(7) 



^ dx 



The dominant term for the even eigenfunctions is 
1 = 4, but we will consider in our calculation even 
terms up to i = 8. From this projection, SZ92 
deduce a fourth-order differential equation for the 
modal function ^j(r), defined in eq (^): 



41)2 ^ Q 



9 d 



pf pvH dr 1 N'^CpT dr 



2d fT d 2o 
dr V P9 dr 



+n^ = o. (8) 



where pi are eigenvalues related to the eigenfunctions 
Fi (for i = 4, p4 = 4.933), x = IGcrT^/^Kp is the ra- 
diative conductivity, A^ the Brunt- Vaisala frequency, 
ri the mean rotation rate of the radiation zone, and 
all the other variables have their usual meaning. 

The problem is simplified by treating the tachocline 
as a boundary layer in which the rapidly varying 
quantity is Qi and all the other physical quantities 
are assumed constant. By doing so, an analytical so- 
lution of the simplified fourth-order differential equa- 
tion is (for simplicity we only give the expression for 
the modal index i = 4): 

QiiC) = fi4(0)\/2exp(-C) cos(C - vr/4) (9) 

where r2j(0) = Qi^, Qi are numerical coefficients de- 
pending on the latitudinal variation of the differen- 
tial rotation and we have defined a non dimensional 
depth: 

C = P^{ncz-r)/d (10) 



with. 



d = r,,,{2^/Nfl\4K/vHfl* 



(11) 



being a parameter related to the tachocline thickness 
h, rbcz the radius and K = x/p^p the radiative diffu- 
sivity, at the base of the convective zone. 

With the latitudinal dependence of the angular ve- 
locity at the base of the convection zone borrowed 
from Thompson et al. (1996): Qbcz/^ir = 456 — 
72x2 _ 422;^ nHz, we have reestimated the coefficient 
Q4 = —1.707 X 10~2^ as well as the ratio between the 
rotation rate of the radiative zone and the equato- 
rial rate: 0,/il.Q = 0.9104. The prediction by Gough 
& Sekii (1997), who consider instead the magnetic 



stresses, is ~ 0.96; presently, the seismic observations 
suggest a rotational ratio of 0.94 it 0.01 (Corbard et 
al. 1999), which is intermediate between these two 
theoretical estimates. 

The inversion of the rotational profile from he- 
lioseismology observations gives the width of the 
tachocline and its location below the convection zone. 
For a comparison between the observational thickness 
and the theoretical one, we have to consider the first 
zero of equation (y), e.g h = 37rd/4/i4 ~ d/2. The 
parameter d is then approximately equal to 2h. In 
the solar models of section y, we will treat h (hence 
d) as an adjustable parameter, in order to agree with 
the helioseismic inversion of the tachocline thickness 
(Basu 1997, Charbonneau et al. 1998, Corbard et al. 
1998, 1999). 

Our purpose now is to obtain an analytical expres- 
sion for the vertical velocity u. We already have an 
expression for the differential rotation, the next step 
is to connect u io Ct. From the conservation of the 
angular momentum, we find a relation between the 
stream function ^ and the differential rotation Vt. 
By writing the streamfunction as ^j = pr'^ui and 
introducing the analytical expression of Cli (eq ^) , we 
finally obtain after integration, the radial dependence 
of Ui; in the case of i = 4 we get: 



the n=4, 6 and 8 terms) 



1 1'Hd 
^ 'bcz 



filQ^exp{-C)cos{C). (12) 



^*4(C) 



From this equation (II4 ), we shall derive an expres- 
sion for the macroscopic diffusivity Dt of the motions 
occurring in the tachocline. 

2.2. Macroscopic diffusion of chemical species 

We determine in this section the influence of the 
mixing on the evolution of chemical species, using 
the expression for the vertical velocity amplitude Ui 
(eq. 1^). The anisotropic diffusion, invoked by SZ92 
to stop the spread of the tachocline deeper in the ra- 
diative zone, interfers with the advective transport 
of chemicals. Chaboyer & Zahn (1992) have shown 
that the result is a diffusive transport in the vertical 
direction, with an effective diffusivity given by: 



Dt 



Uiir) 






(13) 



where Un are the coefficients of the expansion of the 
vertical component of the velocity u in Legendre poly- 
nomials. The eigenfunctions Fi introduced in section 
2.1 may be projected on these Legendre polynomials, 
and to a good approximation one has (keeping only 



Dt 



r 



180 



+ 



4{r) /128y ni(r)' 



546 



V 35 / 1224 
(14) 
One reaches the following expression for the ver- 
tical diffusivity in replacing u^ by equation (|l^ ) in 
equation ([T^, same for n = 6 and n = 8 (c.f Fig ||): 



Dt{C) 



1 1 

1804 



l^H 



(AS)' ..6^2 



\ ^bcz / 



/^4 Q4 



X ex.p{—2() COS (C) + higher order terms. (15) 

where ( is the scaled vertical coordinate defined in 
equation ([To|). 

In this expression, it is assumed that the horizon- 
tal component of the macroscopic diffusivity Dh is 
equal to the horizontal viscosity i^h- This coefficient 
is then introduced in the diffusive part of the equa- 
tion (|l]) for the time evolution of the chemical abun- 
dances Xi, and we could thus obtain solar models 
including a treatment of the tachocline. Fig || com- 
pares the different values of the macroscopic diffusiv- 
ity for models discussed in the next section, to the mi- 
croscopic diffusion (here represented as a horizontal 
solid line). The value of 10 cm?s~^ for Di is deduced 
from the expression of Michaud & Proffitt (1993) for 
^He, in the case of the reference solar model at the 
present age. It must be pointed out that in our calcu- 
lation these microscopic diffusion coefficients as well 
as the microscopic velocity Vi are evaluated for all the 
chemical species present in our gaz mixture, and are 
not assumed constant (see paper I). The amplitude 
of D4Hf, varies only by 10% in the range [0.6,0.72] 
Rq and doesn't change significantly with time. The 
oscillating behaviour of the macroscopic coefficient 
seen in Figure (||) reflects the sub cells structure of 
the tachocline and reduces the lithium depletion in 
comparison with a purely exponential coefficient by 
a factor ~ 1.5-2 and induces stages in the radial chem- 
ical composition profile (c.f 3.2 and FigQ). The mix- 
ing in the tachocline thus depends on the following 
parameters, for a given model of the present Sun: 

• the rotation rate 0,, a measured quantity; 

• the differential rotation rate Ct (or its projection 
Qi) inferred from helioseismology; 

• the Brunt- Vaisala frequency N, determined by 
the extent of convective penetration; 

• the horizontal diffusivity i^h, which presumably 
sides with the differential rotation rate. 

There are at present no firm prescriptions available 
for N and i^h, therefore we must consider these or 
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, n = 8 Bs (- -) and X;L4 -^'^ (t^ii^k -.-). b) Models A (TV = 100 and d = 0.1) (- -), B {N = 25 and d = 0.1) (-.-) 
25 and d — 0.5) (-...-). The horizontal solid line represents approximately the microscopic diffusion Di. 



two combinations of them, as free parameters in our 
calculations. We have chosen to keep N and to take d 
(twice the thickness of the tachocline), because both 
of them are constrained by upper limits drawn from 
helioseismology. 

3. SOLAR MODELS 
3.1. Description of the models 

In paper I (Brun, Turck-Chieze &: Morel 1998), we 
have calculated a new standard model using the CE- 
SAM code (Morel 1997) and have discussed the valid- 
ity of such a model in the light of the present helioseis- 
mic observations. We just recall here the main phys- 
ical inputs present in this reference model. We used 
the solar chemical composition of Grevesse, Noels &: 
Sauval (1996), the OPAL96 opacity tables Iglesias & 
Rogers (1996) and the equation of state of Rogers, 
Swenson Sz Iglesias (1996), the recent updated nu- 
clear reaction rates of Adelberger et al. (1998), the 
prescription of Mitler (1977) for the treatment of the 
nuclear screening reaction rates for a charge Z>2 (Dz- 
itko et al., 1995) and the microscopic diffusion coef- 
ficients suggested by Michaud & Proffitt (1993). 

In the present paper (called paper II), we use 
the "pure microscopic diffusion" model as the ref- 
erence model. We then compute solar models using 
the macroscopic coefficient Dt described by equation 
(p^) and follow the time evolution of these solar mod- 
els from the pre-main sequence (PMS) until an age 
of 4.55 Gyr, making the age of our present Sun 4.6 
Gyr (including about 0.05 Gyr of pre-main sequence). 
The adaptable time step increases from ~20 years up 
to 200 Myr during the evolution. We start from the 



PMS so as to properly follow the ^Li abundance evo- 
lution, including the macroscopic diffusion coefficient 
along all the solar phases. We introduced an updated 
nuclear reaction rate for ^Li(p,a)^He, Engstler et al. 
(1992) instead of the Caughlan & Fowler (1988) one. 
The results of the solar models we calculated for this 
study are summarized in Tables || and ^ and the pho- 
tospheric abundances of the elements which are sensi- 
tive to the adding macroscopic process are compared 
to observations in Table ^. 

3.2. Time independent rotation 

We first compute solar models including the effec- 
tive diffusivity Dt (eq. |I5|), with a tachocline thick- 
ness h of 0.05 Rq {d = 0.1 Rq) in agreement with 
the helioseismic measurements (Basu 1997, Corbard 
et al. 1998, 1999, Charbonneau et al. 1998), a Brunt- 
Vaissala frequency A^ of 100 or 25 /uHz to take into 
account the rapid variation of this quantity at the 
base of the convective zone, and a rotation fi at the 
base of the convective zone of 0.415 //Hz; we call these 
coefficients A and B respectively (see Table |l]). The 
reference model has a present photospheric helium 
abundance of 0.2427 in mass, corresponding to an 
^He diffusion of 10.8%. This value of Y^ is slightly too 
low if we compare with the Basu & Antia (1995) value 
for the OPAL equation of state, Y^ = 0.249 ± 0.003. 
When introducing the macroscopic diffusion coeffi- 
cient, we mix back helium in the convection zone, 
inhibiting the microscopic diffusion by 15% or 25% 
and leading to a present photospheric ^He abundance 
Ys = 0.2452 or Y^ = 0.2473 respectively for models 
A and B (c.f Tables |^ and y, see also Brun, Turck- 
Chieze & Zahn 1998). This is in better agreement 



6 



with the hehoseismic values. The decrease of the 
Brunt- Vaisala frequency N leads to an increase of 
the efficiency of the mixing in the tachocline layer 
and reduces the settling of the chemical species. We 
also note that the initial helium content Yq is smaller 
than 0.27 in mass (cf Table || and also paper I). 

In Figure ^ we present the difference between the 
initial helium composition Yq and the present photo- 
spheric one Yg for several models including the ref- 
erence one. Before the introduction of microscopic 
diffusion, this difference was null because the present 
photospheric values were adopted as the initial com- 
position of the star. However now, when we compute 
a standard model including the microscopic diffusion 
and the settling of chemical species, the initial abun- 
dances are assumed to be different from their present 
day values and there is an iteration on the composi- 
tion to get the photospheric Z/X ratio in agreement 
with the one measured at the present time: Z/X= 
0.0245 ± 0.003. 

The prominent feature of Figure || is the smoother 
composition profile for models including the macro- 
scopic diffusivity Dt, in comparison with the refer- 
ence model, which presents a sharp composition gra- 
dient just below the convective zone. The compo- 
sition profile is also almost constant below the con- 
vective zone over a distance which is related with the 
value of d and A'^. A smaller Brunt- Vaissala frequency 
value induces a greater mixing, and then an extended 
plateau. The characteristic steps in the radial pro- 
file come from the cosine term in the expression of 
Dt (eq. [l|) (see also EUiott (1997)). An explana- 
tion could be that the simplified assumption in the 
treatment of the tachocline of SZ92, of a constant 
turbulent horizontal viscosity vh, induces this radial 
dependence of the coefficient Dt and thus the steps in 
the composition profile. A more careful study should 
be done to verify this statement. 
Concerning the Li and ^Be depletion, we see in Ta- 
ble ^ that the mixing increases the lithium depletion 
by a factor up to ~ 4, without burning ^Be, which 
agrees with the revised determination by Balachan- 
dran & Bell (1998). Nevertheless, such an increase 
of the ^Li depletion is still insufficient to reproduce 
the photospheric lithium depletion (Grevesse, Noels 
& Sauval 1996, Cayrel 1998 and references therein). 

3.3. Time dependent rotation 

In this section we shall take into account the time 
evolution of the tachocline width and the efficiency 
of the mixing from the early phase of the Sun until 
now. In the previous section, we have used a constant 
value for d based on the present observations of the 
tachocline thickness in the effective diffusion coeffi- 
cient Dt- Thus, in order to improve the description 



of tachocline mixing, we now introduce a time de- 
pendence of the coefficient Dt and of the tachocline 
thickness h (hence d), related to the global and the 
differential rotation of the Sun as follows: 



Dt oc I'H 



_d_ 



Q\ oc ^v^^iSll^f , d oc fi^/V^^l/^ 

(16) 
where we have used equations ([T5| ) and (11). As- 
suming that the turbulent viscosity is proportional 
to the differential rotation (i.e vn oc Q) and intro- 
ducing the dependence of differential rotation on ro- 
tation observed by Donahue, Saar & Baliunas (1996) 
(A oc f70.7±o.i)^ ^g £j^j^iiy Q^x.wi for Dt{S^) and d(0) 
the following scaling laws: 

Dt oc f70.75±0.25 ^ d OC f](l-3T0-l)/l (17) 

We conclude that the tachocline mixing was stronger 
in the past both because that layer was thicker and 
because the diffusivity was larger. Note that the ob- 
servational uncertainties in the relation between 
and J7 act in opposite way on the dependence of 
Dt{^) and d{i^) (i.e a stronger mixing implies a 
smaller tachocline extent). We render the mixing 
in the tachocline time-dependent, through Z?j'(r2(t)) 
and d{^{t)), by using the spin-down law Q{t) oc t^^''^ 
which was deduced by Skumanich (1972) from the 
rotation rate of stellar clusters of different ages. 

Models At, Bt and Ct correspond to solar evolu- 
tions with time-dependent coefficients given by equa- 
tion (ll^) and different values for A^ and the present 
tachocline thickness d (see Table 1 and 2). The first 
two cases correspond to the same values of d and A^ 
that existed in the previous section but with a time 
dependence of the coefficient Dt- The third shows 
the influence of a thinner tachocline {d = O.OSi?©) 
as suggested by more recent studies (Corbard et al. 
1999). The introduction of this time-dependence in- 
creases the inhibition of the microscopic diffusion of 
helium during the evolution in comparison with the 
previous solar models considered in section ^T^ (see 
Table ^ . These models show an inhibition of ^He up 
to 27% for model Bt {d = Q.l, N = 25) leading to a 
surface helium Y^ = 0.2477, very close to the seismic 
value of Basu & Antia (1995) . 

We plot the ^Li and ^Be radial profiles in Figure 
1^, for the reference model and the time dependent 
models calibrated in Z/X at 4.6 Gyr, with the corre- 
sponding radial temperature scale. These two light 
chemical elements are sensitive to mixing processes 
occuring in the stars because their nuclear burning 
temperatures are very low (resp. 2.5 10^ and 3.2 10^ 
K) (Baglin & Lebreton 1990). It is well known that 
the solar ^Li is depleted by a factor ~ 100 (Grevesse, 
Noels & Sauval 1996) and recent revised abundance 
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Fig. 2. — Radial profile of the difference of He composition between the initial and present values for the reference solar model 
including only microscopic diffusion (solid line) and solar models where we add a macroscopic mixing due to the presence of the 
tachocline: coefficient A (- -) and B (-.-) (see section B.3 and Tables hj, H for the model characteristics). 



determinations have shown that ^Be abundance is 
practically not depleted (Balachandran & Bell 1998) 
contrary to what was thought previously. These new 
observational constraints can only be satisfied if those 
chemical species are mixed in a rather thin layer be- 
low the convective zone, which is the case for the 
tachocline mixing considered here. 

We clearly see in Figure ^ that this mixing process 
modifies the distribution of lithium but not of beryl- 
lium (except the flat plateau for the mixed models in 
comparison with the "pure" diffusion one) . These re- 
sults are confirmed in Table S where we show the 
initial over present ratio of the ^Li and ^Be. We 
also remark that models including a macroscopic dif- 
fusivity with a low value of the Brunt- Vaisala fre- 
quency A^ = 25 (i.e Bt and Ct) burn more ^Li, be- 
cause of their higher absolute value (> 10^ cm^s~^ 
comparison with the microscopic diffusion one ( 
cm^s~^), than models with a higher N = 100 (i.e At) 
(Figure 1). The models Bt and Ct strongly inhibit 
the settling of the chemical species and mix all the 
elements over a distance where the lithium burning 
temperatures are reached r ~ O.67i?0. The charac- 
teristic time for this mixing, for a Dt of 10^ cm^s~^ 
is ~ 3 Myr. This time is intermediate between the 
convective cell turnover time (~ 1 month) and the 
microscopic diffusion one (> 10 Gyr). We also see in 
Table |3| that the time dependent diffusivity improves 
the ^Li depletion, where a value of ~ 100 is obtained 
without destroying ^Be and without increasing the 
^He/^He surface ratio over the past 3 Gyr, as de- 
duced by Geiss & Gloeckler (1998) from meteorites 
and solar wind abundance measurements. As already 
mentioned these abundance constraints can only be 
satisfied if the mixing occurs in a thin layer below the 



m 
~ 10 



convective zone, which is the case for the tachocline 
considered here. 

In Figure |^ we present the lithium depletion occur- 
ring during the Sun evolution for the different mod- 
els presented, plus one without microscopic diffusion. 
Clearly, only the diffusion models including mixing 
in the tachocline show a substantial depletion during 
the main sequence evolution as it is observed in dif- 
ferent clusters. One notes that the lithium depletion 
is largely dependent on the Brunt Vaissala frequency 
adopted at the base of the convective zone and on 
the time dependence of the macroscopic coefficient. 
The modification of the tachocline thickness has only 
little impact on the lithium depletion as the higher 
value of the coefficient largely compensates the re- 
duction of the thickness. However one sees that the 
^Li depletion slopes in these two models (i.e Bt and 
Ct) are not the same. The strong time dependent 
mixing with A^ = 25 presents a correct value of the 
solar Li depletion. The lithium depletion during the 
PMS is probably overestimated due to the crude spin- 
down law we have adopted. A more detailed analysis 
of these phases is under study (see also Vauclair & 
Richard 1998 and Montalban & Schatzman 1996). 



4. DISCUSSION 



4.1. 



The sound speed profile 

We now compare the sound speed profile of the 
different models with the helioseismic inversion re- 
sults obtained with the GOLF-I-MDI data (Lazrek 
et al. 1997, Roca Cortes et al. 1998, Kosovishev 
et al. 1997) and described in Turck-Chieze et al. 
(1997, 1998) (Figl), we see that the peak localized 
at r ~ 0.68Rq is practically erased by this process, 
depending on the strength of the mixing. 



T6(K) 




r/Ka 



Fig. 3. — ^Li and ^Be radial profile for several models (superimposed the corresponding temperature Tg =T/10^ K obtained 
for the reference model). Respectively: reference (solid line), models At (dash), Bt (dash dot) and Ct (dash three dots). The 
mixing below the convective zone modifies only the ^Li profile and not the ^Be one (except the flat plateau for the mixed models). 
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Fig. 4. — Time dependent depletion of ^Li for several solar models: no microscopic diffusion ( ), with microscopic diffusion 

(solid line) and with mixing in the tachocline thickness: models A (-.-), B (-...-), then with time dependent mixing, models At (...) 
and Bt (long dash), Ct (thick -.-) and Btz (Zq =Zo*'* = 0.01959) (thick - -). We superimposed on the theoretical curves the open 
cluster observations: {aPer: Balachandran et al. (1996); Pleiades: Soderblom et al. (1993b); UMaG: Soderblom et al. (1993a); 
Hyades: Thorburn et al. (1993); NGC752: Balachandran (1995). (Figure adapted from Vauclair & Richard 1998 and cluster age 
uncertainties deduced from Lebreton et al. 1997 and private communication.) 



In all the models shown in this figure (El), we have 
adjusted the initial composition in order to get at the 
solar age, the present observed photospheric heavy 
element/hydrogen ratio, i.e Z/X=0.0245. This cali- 
bration causes the swell of the small bump located 
between the nuclear core and the convective zone, 
since it reduces the initial metallicity (see Table 2). 
It is well known that the heavy element abundance 
has a direct infiuence on the solar structure, mainly 
through the opacity coefficients and hence it influ- 
ences the sound speed profile between 0.4 and 0.7 Rq 
(Turck-Chieze 1995, 1998 and paper I). The result 
is quite different if we do not constrain the photo- 



spheric Z/X value of the mixed model, starting with 
the heavy element composition Zq = 0.01959 of the 
reference model for example, we obtain a different 
chemical composition at the solar age (see for exam- 
ple columns 3 and 4 of Table y for the coefficients Bt 
and Btz)- In Figure |^, we draw the corresponding 
squared radial sound speed differences between the 
models and the induced solar value for two models 
where we have kept (i.e coefficient Bt) or relax {Btz) 
this surface heavy element constraint. We see that in 
both cases where the macroscopic term is introduced, 
the bump around 0.68 Rq is practically erased by the 
introduction of the tachocline mixing, but the cali- 
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Fig. 5. — Sound speed square difference between GOLF+MDI data, the reference model (solid line) and models including a 
macroscopic term: At ( ), Bt (-.-), Ct (-..•-) (see parameters in Table 2). 



brated model is modified along the whole structure 
and not only close to the bump. It is not the case 
for the non-calibrated model where Zq is fixed. This 
is partly illustrated by Figure |^, where the non cal- 
ibrated model changes only locally the structure of 
the solar model (dash dot line). 

A similar result is obtained by Elliott, Gough & 
Sekii (1998) invoking an advection of the heavy ele- 
ments at the base of the convection zone. They de- 
duced from the sound speed profile a low value of the 
tachocline thickness (h = 0.02) which could be due 
to the macroscopic coefficient used. In this estimate, 
they did not calibrate their model in Z/X and ob- 
tained a result similar to ours (dash dot line in figure 

W)- 

4.2. Solar composition 

Examining the different observables, it seems that 
the model using the coefficient Btz (see Tables § and 
HI) corresponds today to a good representation of the 
present Sun. The agreement on the sound speed is 
significantly improved and the observed photospheric 
light elements (^He, ^He, ^Li, ^Be) are well repro- 
duced. In this model, two effects contribute to this 
improvement; the introduction of the mixing at the 
base of the convection zone, which partly inhibits the 
microscopic diffusion, and the induced increase of the 
heavy elements at the surface and consequently in 
the intermediate region due to this inhibition. This 
model treats the tachocline instability as a local phe- 
nomenon. 

In Figure ^b we show the heavy element profile 
along the solar radius of the reference model and 
the mixed models (coefficients Bt & Btz) with and 
without the Z/X=0.0245 calibration. We clearly see 



that without the Z/X calibration there is a difference 
in the Z proffie, corresponding to an increase of the 
opacity and then to a reduction of the disagreement 
between the seismic sound speed and the theoretical 
one. 

In paper I, we had already mentioned that the 
sound speed difference between the Sun and the 
model may suggest that a slight increase of the heavy 
element composition of about 5% is required. One 
notes that the results we obtain without specific cali- 
bration but keeping the initial heavy element com- 
position of the reference model, lead to a photo- 
spheric heavy element composition of Z/X=0.0255. 
This is still consistent with the present observa- 
tional uncertainties (about 10%) even if a very re- 
cent work seems to suggest a possible reduction of 
the CNO abundances rather than an slight increase 
(Grevesse & Sauval 1998). We have verified this 
statement by performing a "pure" microscopic model 
with Z/X=0.0255. If one compares the sound speed 
proffie of model Btz (dot dash line) to the reference 
model with Z/X=0.0245 (Fig |6|), it shows a better 
overall agreement with the seismic data. Otherwise, 
if one compares this model to a "pure" microscopic 
model with Z/X=0.0255, one reaches a different con- 
clusion for the radiative part (identical to the com- 
parison of model Bt with the reference model Fig ^). 
This intricate process shows that it is not so easy to 
find the real behaviour of the Sun. If we increase 
the initial heavy element abundance a little more, we 
destroy the agreement of the tachocline region and 
obtain a photospheric composition which is not in as 
good agreement with the present observations. So 
extensive studies are surely useful to properly disen- 
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Fig. 6. — Sound square speed difference between GOLF-fMDI data, the reference model (solid line) and two models including 
a macroscopic term: Bt calibrated in Z/X (- -) and Btz with a non calibrated Zo —TJi}'^ — 0.01959 (-•-). 



tangle the competing processes. 

The remaining part of the sound speed discrepancy 
located between the nuclear core and the convection 
zone (see Fig g) may come from uncertainties in the 
equation of state or opacity coefficients (Turck-Chieze 
et al. 1998), or from our simplified treatment of the 
microscopic diffusion. In the present study, we do 
not consider the different states of ionisation of the 
different species or the radiative acceleration in the 
calculation of the microscopic diffusivity (Turcotte et 
al. 1998). We also cannot exclude the presence of 
other mixing processes, for example, the meridional 
circulation or the effect of the internal waves or mag- 
netic field. It seems difficult to think that the residual 
discrepancy comes from insufficient knowledge of the 
element abundances, as an increase of the heavy ele- 
ment composition will have a dramatic effect in the 
region of the energy transport transition. 

4.3. Nuclear reaction rates 

The depletion of lithium is of course dependent on 
the nuclear reaction rates used. In this study, we 
have used the recent coinpilation of the NACRE con- 
sortium (Angulo et al. 1999), which recommends the 
nuclear reaction rate ^Li(p,a)'^He of Engstler et al. 
(1992). The astrophysical factor is increased by 30% 
in comparison with the older value recommended by 
Caughlan & Fowler (1988), mainly due to a higher 
accuracy numerical integration. The updated value 
of 5*71(0) = 0.0593 Mev barn increases the lithium de- 
pletion during the Sun evolution by ~ 40% (see Table 
3 model Bto and Bt). This effect is introduced in all 
the results presented in this study and constitutes 
an important point which must be considered in the 



lithium depletion problem in stars. It is interesting 
also to notice that the reaction rate of the ^Li(p,7)^B 
which is not included in our study, mainly because 
this cross section was considered as very small, has 
been reestimated these last years and found 100 times 
greater than in Caughlan & Fowler (1988). The im- 
pact on the present study must be smaller than 10%, 
considering the relative importance of the two reac- 
tions. 

4.4. Neutrino predictions 

Of course, one issue of improving solar models is to 
predict more correctly the neutrino fluxes produced 
in the central region of the Sun (Turck-Chieze 1999). 
Tables 1 and 2 show the predicted values for all the 
models considered. One notices a very small depen- 
dence on the macroscopic coefficient used, but nev- 
ertheless two kinds of results depending on whether 
the models are calibrated in Z/X or not. This dif- 
ference is connected to the effect of the inhibition 
of the microscopic processes which leads to a slight 
decrease of the neutrino fluxes. Indeed, we observe, 
in the calibrated models, a reduction in comparison 
with the reference model of -1.6% for the Gallex and 
Sage experiments (resp. Hampel et al. 1996 and Ab- 
durashitov et al. 1996), of -5.1% for the Homestake 
experiment (Cleveland et al. 1998) and -5.6% for Su- 
perKamiokande (Suzuki et al. 1998). It is useful to 
recall that the general effect of the introduction of the 
microscopic diffusion in solar model was an increase 
of about 20% on the ^B neutrino flux prediction, so 
this corresponds to 25% of this effect. But, one may 
also notice that the effect of inhibition could be ex- 
actly compensated by the effect of the increase of the 
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a macroscopic term: Bt calibrated in Z/X ( ) and Btz with a non calibrated Zo =Zo* = 0.01959 (dash dot line), b) Heavy 

Element Composition for the reference model including only microscopic diffusion (solid line) and two models with macroscopic 
mixing due to the tachocline: Bt calibrated in Z/X 
the model characteristics). 



and Btz with a non calibrated Zo =Zg" = 0.01959 (-.-). (see Table | for 



heavy elements by 5% (reference model and model 

These effects are too small to solve the neu- 
trino problem. This is hardly a surprise, since the 
tachocline mixing occurs at the base of the convec- 
tion zone, hence too remote to have a serious impact 
on the structure of the nuclear core. Even if this 
study is not crucial for the neutrino fluxes it gives 
certainly a proper description of the chemical trans- 
port of elements in the external part of the Sun and 
definitively excludes a large increase of the neutrino 
fluxes by this kind of processes as was thought five 
years ago. The introduction of a macroscopic diffu- 
sivity also destroys the criticisms on the impact of 
the pure microscopic diffusion which was well known 
to be partly inhibited by some turbulent processes 
and which seems now well under control, at least as 
far as the extension in the Sun is concerned. We be- 
gin to include dynamical effects properly in stellar 



models and one cannot exclude that the cumulative 
effect of different modifications of the physical inputs 
(microscopic and macroscopic) may lead to smaller 
neutrino fiuxes. We will not comment, in this paper, 
on the discrepancy of the sound speed in the cen- 
tral region because this is more difficult to estimate 
(Turck-Chieze & Basu 1998, Gonczi et al. 1998) and 
to interpret there, due to the sensitivity of the sound 
speed to a lot of competing processes (see paper I). 
Works are in progress to have good constraints on the 
core. 

5. CONCLUSION 

In this paper, we have shown that the presence of 
mixing in the tachocline layer partly inhibits the mi- 
croscopic diffusion process and improves the agree- 
ment with the helioseismic and photospheric abun- 
dance data. Our results make use of updated he- 
lioseismic constraints on the present extent of this 
layer, and they are based on a physical treatment of 
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the tachocline layer (SZ92) which involves only one 
additional parameter, namely the Brunt- Vaisala fre- 
quency (which we assume constant in our model). 
The solar models thus obtained are quantitatively 
in good agreement with the observational data and 
our favoured model is model Btz (cf. Tables 2 &: 
3). When introducing a time dependent tachocline 
mixing, it is possible to reach the observed Li pho- 
tospheric abundance at the present solar age, without 
destroying ^Be or bringing too much ^He at the sur- 
face. This does not necessarily mean that our mixing 
process is the only one which operates in the Sun: 
others may contribute, such as mixing by gravity 
waves (Montalban & Schatzman 1996) or meridional 
circulation (Zahn 1992, Vauclair & Richard 1998). 
But the lack of depletion of ^Be proves that the mix- 
ing occurs in a rather thin layer below the convective 
zone, which is the case in our models. The model 
which seems the most promising assumes a Brunt- 
Vaisala frequency A^ = 25/xHz; it leads to a slight 
increase of the heavy element abundance by 5% at 
the surface, compatible with the present uncertain- 
ties. A detailed study of the feedback of such mixing 
on the structure of the Sun reveals that its effect is 
not confined to the base of the convection zone, but 
that it spreads into the region below, except if one 
relaxes the fixed Z/X constraint and let the model 
adjusts within the heavy element abundance incerti- 
tudes (model Btz)- 



Prom the results of papers I and II, we could deduce 
a galactic enrichment of ^He. We find that the so- 
lar ^He primordial abundance is Yq = 0.27]'^o!ooi ^^"^ 
leads to an enrichment of 10.6Ti'3% ^oi^ the value 
of 0.244 ± 0.002 deduced by Thuan & Izotov (1998) 
from observations of a large sample of low-metallicity 
H II regions in blue compact dwarf galaxies. 

This first attempt to calibrate dynamical processes 
with helioseismology data from SOHO leads to the 
idea that the neutrino problem could be seen from 
a new angle, by considering reliable solar models 
with macroscopic processes, never implemented pre- 
viously. The process invoked here concerns the 5% 
part of the external mass and consequently has little 
impact on neutrino predictions. But, even if the solar 
models are close to the seismic data for r < 0.2Rq, 
the remaining discrepancies are still significant (~ few 
a) in comparison with seismic uncertainties, and they 
suggest that other physical processes, neglected so 
far, may operate there. 

To conclude, this study encourages the introduc- 
tion of new processes in stellar structure theory, be- 
yond the standard model framework. Of course we 
are aware of the crudeness of the present model: a 
more refined treatment would include the variation 
with depth of the horizontal turbulent diffusivity and 
of the Brunt- Vaisala frequency, as well as a more pre- 
cise spin-down law, and it would not ignore the mag- 
netic field. 
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Table captions 



Definition of the parameters used in Tables || and |2|: 

d: is twice the tachocline thickness h, N: Brunt- Vaisala frequency, a: mixing length parameter, Yq, Zq, 
(Z/X)o: initial helium, heavy element and ratio heavy element on hydrogen, Y^, Z^, (Z/X)^: idem for surface 
compositions, Hbzc, T^bcz are the radius and temperature at the base of the convective zone, Yc, Zc, Tc, pc'- 
central helium, heavy element contents, central temperature and density, Ga, ^^Cl, ^B respective neutrino 
predictions for the gallium, the chlorine and water detectors. OPAL/ A means that we use OPAL96 opacities 
and Alexander (1994) for low temperature. 



Table 1 
Solar models with a time independent macroscopic diffusivity describing the tachocline mixing and 

COMPARED TO THE REFERENCE MODEL, EACH OF THEM CORRESPONDS TO DIFFERENT VALUES OF THE BRUNT 

VAISSALA frequency and to the TACHOCLINE THICKNESS. 



Parameters 


Ref 


A 


B 


Opacities 


OPAL/A 


OPAL/A 


OPAL/A 


Diffusion 


yes 


yes 


yes 


Age (Gyr) 


4.6 


4.6 


4.6 


d {t/Rq) 


- 


0.1 


0.1 


N (/uHz) 


- 


100 


25 


(Z/X). 


fixed 


fixed 


fixed 


a 


1.766 


1.753 


1.745 


Yo 


0.2722 


0.2703 


0.2691 


Zo 


0.01959 


0.01919 


0.01893 


(Z/X)o 


0.0277 


0.0270 


0.0266 


Y, 


0.2427 


0.2452 


0.2473 


Is 


0.01803 


0.01797 


0.01792 


(Z/X)« 


0.0245 


0.0245 


0.0245 


Rfczc /R-o 


0.7133 


0.7145 


0.7154 


Tb.e X 106 (K) 


2.190 


2.181 


2.174 


Y, 


0.6405 


0.6382 


0.6367 


Zc 


0.02094 


0.02051 


0.02024 


Tc X 10^ (K) 


15.71 


15.68 


15.66 


Pc [g/cTD?) 


153.1 


152.9 


152.7 


^^Ga (SNU) 


127.1 


126.0 


125.2 


3^01 (SNU) 


7.04 


6.84 


6.71 


8B (lOVcmVs) 


4.99 


4.83 


4.73 
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Table 2 
Solar models with a time dependent macroscopic diffusivity describing the tachocline mixing and 

COMPARED TO THE REFERENCE MODEL, EACH OF THEM CORRESPONDS TO DIFFERENT VALUES OF THE BRUNT 
VAISSALA frequency and to the tachocline thickness, the model Btz SHOWS THE EFFECT OF NON 

CALIBRATED Z/X (SEE TEXT). 



Parameters 


Rcf 


At 


Bt 


Bt. 


Ct 


Opacities 


OPAL/A 


OPAL 


OPAL/A 


OPAL/A 


OPAL/A 


Diffusion 


yes 


yes 


yes 


yes 


yes 


Age (Gyr) 


4.6 


4.6 


4.6 


4.6 


4.6 


d (r/Re) 


- 


0.1 


0.1 


0.1 


0.05 


N (/zHz) 


- 


100 


25 


25 


25 


(Z/X), 


fixed 


fixed 


fixed 


free 


fixed 


a 


1.766 


1.752 


1.743 


1.755 


1.748 


Yo 


0.2722 


0.2701 


0.2689 


0.2722 


0.2695 


Zo 


0.01959 


0.01914 


0.01889 


0.01959 


0.01903 


(Z/X)o 


0.0277 


0.0269 


0.0265 


0.0277 


0.0267 


Y. 


0.2427 


0.2455 


0.2477 


0.2510 


0.2464 


z. 


0.01803 


0.01796 


0.01791 


0.01858 


0.01794 


(Z/X), 


0.0245 


0.0245 


0.0245 


0.0255 


0.0245 


R.bzc/R.0 


0.7133 


0.7146 


0.7155 


0.7141 


0.715 


Tfc,, X 106 (K) 


2.190 


2.180 


2.173 


2.194 


2.177 


Yc 


0.6405 


0.6379 


0.6365 


0.6405 


0.6372 


Zc 


0.02094 


0.02047 


0.02019 


0.02094 


0.02034 


Te X 10^ (K) 


15.71 


15.67 


15.65 


15.71 


15.66 


Pc (g/cm^) 


153.1 


152.8 


152.7 


153.1 


152.8 


^iQa (SNU) 


127.1 


125.9 


125.1 


127.1 


125.5 


37C1 (SNU) 


7.04 


6.82 


6.69 


7.04 


6.76 


8B (lOVcmVs) 


4.99 


4.81 


4.70 


4.99 


4.76 



Table 3 
Surface abundance variation of ^He/**He during the last 3 Gyr, surface abundances of ^He and 

HEAVY elements Z, AND ABUNDANCE RATIO INITIAl/SURFACE FOR ^Ll AND ^Be FROM OBSERVATIONS AND FOR 
SOLAR MODELS AT THE SOLAR AGE (SUBSCRIPTS t FOR TIME DEPENDENT MODELS, tz FOR TIME DEPENDENT MODEL 
WITH Zo —T^Q^^ ~ 0.01959 AND to FOR TIME DEPENDENT MODEL WITH CAUGHLAN & FOWLER (1988) '^Li(p,Q!)'*HE 

NUCLEAR REACTION RATE) . 





Obs 


Rcf 


A 


B 


At 


Bt 


Bt, 


Bto 


Ct 


d {r/R-Q) 


<0.1 


- 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.05 


N (mHz) 


- 


- 


100 


25 


100 


25 


25 


25 


25 


(3He/4Hc), 


max 10% 


2.28% 


2.14% 


2.01% 


2.11% 


2.0% 


2.02% 


2.0% 


2.07% 


■^He, 


0.249±0.003 


0.2427 


0.2452 


0.2473 


0.2455 


0.2477 


0.2509 


0.2477 


0.2464 


(Z/X). 


0.0245±0.002 


0.0245 


0.0245 


0.0245 


0.0245 


0.0245 


0.0255 


0.0245 


0.0245 


^Lio/^Li, 


-100 


-6 


-8 


-22 


-12 


-91 


-134 


-64 


-89 


^Bco/'^Bc, 


1.10±0.03 


1.115 


1.093 


1.086 


1.093 


1.118 


1.125 


1.118 


1.092 



